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Abstract 


The  final  results  and  academic  achievements  are  reported  on  a program 
aimed  at  investigation  of  external  burning  propulsion.  Concentration  was 
aimed  at  bluff-base,  axisymmetric  bodies  in  supersonic  flight  with  an 
approach  boundary  layer  that  is  fully  turbulent.  An  adiabatic  near  wake 
theory  was  developed  which  is  in  excellent  agreement  with  experiment. 

A parametric  investigation  of  external  burning  performance  showed  that 
while  high  thrust  levels  can  be  achieved  the  specific  impulse  is  lower 
than  expected.  Consequently,  future  studies  should  concentrate  upon  a 
combination  of  base  burning  and  external,  burning. 
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Introduction 


There  arc  a great  number  of  potential,  air-to-air,  air-to-ground  and 
ground-to-air  wcaponi'y  missions  that:  require  either  a sustain  or  mild 
acceieration  phase  of  the  missile  trajectory  or  which  could  benefit  from 
a substantial  drag  reduction  during  a portion  of  the  trajectory.  Further- 
more, many  of  these  missions  require  operation  sufficiently  low  in  the 
atmosphere  that  airbreathing  propulsion,  if  it  is  competitive  with  the 
rocket  is  attractive. 
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Investigated  in  this  program  was  the  concept  of  external  burning.  In 
its  purest  form  this  consists  of  burning  outside  of  all  shear  layers  near 
the  body,  in  the  adjacent  supersonic  stream.  Moreover,  this  program  concen- 
trated upon  external  burning  behind  axisymmetric  bodies  with  a turbulent 
approach  boundary  layer.  The  idea  is  that  compression  waves  focusing  upon 
the  near  wake  will  be  able  to  transmit  the  high  pressure  through  the  sub- 
sonic portion  of  the  wake  and,  hence,  raise  the  base  pressure.  This  concept 
is  not  to  be  confused  with  base  burning,  whereby  a combustible  is  entrained 
into  the  near  wake.  While  base  burning  is  also  promising  for  base  pressure 


elevation,  the  concentration  on  external  burning  was  motivated  by  the  desire 
for  high  base  thrust  levels. 

The  program  was  analytical  in  nature  but  used  data  from  other  programs  to 
validate  the  theory  as  it  was  developed.  The  program  proceeded  in  two  steps. 
The  first  was  to  create  an  accurate,  but  computationally  fast,  theory  for 
the  near  wake.  The  second  was  to  incorporate  external  burning  into  the  theory. 
The  resuLts  of  the  first  effort  are  attached  as  Appendix  A and  of  the  second 
as  Api'indix  B. 
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Summary  of  I'esu  Lis 

An  integral  method  was  used  to  construct  a theory  of  the  near  wake 
behind  bluff  base  bodies  in  supersonic  flight  with  a turbulent  approach 
boundary  layer,  The  method  is  computationally  fast  and  has  been  checked 
against  numerous  experiments,  with  excellent  agreement  in  base  pressure, 
wake  length  scales,  and  detailed  field  quantities,  Iho  following  conclu- 
sions may  be  drawn  from  the  theory: 

(1)  'Iho  solution  of  the  inviscid  supersonic  flow  adjacent  to  the.  near 
wake  plays  a strong  role  in  the  overall,  solution.  It  must  be  treated 
in  a nearly  exact  (method  of  characteristics)  manner. 

(2)  The  base  pressure  is  a strong  function  of  flight  Mach  number  and 
decreases  with  a Mach  number  increase.  However,  detailed  flow  field 
quantities  such  as  the  location  of  the  rear  stagnation  point  and 
velocity  on  the  dividing  streamline  are  only  weak  functions  of  Mach 
number, 

(3)  The  base  pressure  is  a weak  function  of  upstream  boundary  layer 
thickness.  However,  the  detailed  fLow  field  quantities  are  strongly 
affected. 


(4)  One  area  where  there  is  poor  agreement  between  theory  and  experiment 
is  that  of  base  bleed.  The  theory  underpredicts  the  value  of  base 
pressure  for  a given  bleed  rate. 

(5)  The  integral  method,  which  used  the  boundary  layer  approximation 
for  the  viscous  near  wake,  ran  into  some  ’singularity1  problems  with 
some  velocity  profiles  and  governing  equations.  More  investigation  is 
warranted  to  determine  the  nature  of  these  singularities,  and  to  sys- 
tematize the  procedure  for  overcoming  these.  Experiments  show  that, 


mummim 
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in  fact,  the  boundary  layer  approximation  is  i>ooi'  near  the  base  plane; 
this  assumption  should  be  relaxed  in  future  treatment  of  tile  problem. 

The  above  wake  theory  was  then  combined  with  a theory  incorporating 
combustion  external  to  the  near  wake.  The  combustion  zone  was  treated  as 
a nearly  one -dimensional  premixed  region  with  a centrifugal  correction. 

The  following  conclusions  may  bo  drawn. 

(1)  'flic  base  pressure  may  be  elevated  to  a net  thrust  condition  by 
this  method,  but  the  specific,  impulse  performance  with  state-of- 
the-art  fuel  rich  propellants  is  poor;  it  may  be  beaten  by  a 
conventional  rocket.  Consequently,  it  is  believed  that  the  pure 
external  burning  concept  should  not  bo  considered  aLone,  but  in 
conjunction  with  base  burning.  External  burning  provides  the  means 
for  high  thrust,  whereas  base  burning  can  provide  high  efficiency 
at  lower  ba.se  pressure  rises. 

(2)  The  primary  variables  that  affect  external  burning  performance  are 
Mach  number,  fuel-air  ratio,  total  combustible  mass  flow  and  the 
fuel  calorific  value,  flic  performance  is  insensitive,  within  reason- 
able tolerance,  to  heating  zone  placement  and  heat  release  distri- 
bution* 


(3)  It  is  imperative  that  a theory  be  constructed  for  a combined  base 
burning  and  external  burning  system  in  order  that  an  upper  limit 
of  performance  can  be  established* 
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tig.  2 Centerline  pressure-  variation  with  axial  distance  and  com- 
parison with  experiments  near  Mach  a;  0 = rear  stagnation  point,  C 
= Crocco-I.ees  critical  point,  l>  - wake  neck.  E - centerline  Mach  1 
point. 
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Figures  2 and  3 show  a comparison  of  some  of  the  results  of 
the  present  theory  with  the  experimental  data,  and,  as  seen, 
good  agreement  is  obtained.  However,  theory  seems  to 
slightly  overprcdict  the  base  pressure  at  large  Mach  numbers 
( - above  2.5)  and  slightly  ttndcrpredicl  the  base  pressure  at 
low  Mach  numbers  (-  below  1.8).  It  is  found  that,  if  the  eddy 
viscosity  is  multiplied  by  a factor  x/Af’,/2,  a much  better 
correlation  with  the  experimental  results  can  be  obtained  for  a 
wider  range  of  Mach  number.  Also,  the  variation  of  centerline 
Mach  number  and  the  effect  of  boundary-layer  thickness  on 
base  pressure  compared  well  with  experiments.  These  com- 


parisons with  the  experiments  thus  show  the  adequacy  of  the 
present  modeling  of  the  corner  region  and  the  shear  stress.  A 
number  of  studies,  using  the  present  program,  were  made  and 
the  major  conclusions  arrived  at  areas  follows: 

1)  The  solution  of  ouicr  regime  plays  a very  important  role 
in  the  present  problem:  the  replacement  of  AMC  with  the 
Prandtl-Meyer  relation  resulted  in  very  poor  results. 

2)  For  better  prediction  of  flowfield  details,  a good  velocity 
profile  is  essential.  However,  base  pressure  is  affected  in  a 
minor  way  because  of  a change  in  profiles. 

3)  The  base  pressure  is  a strong  function  of  upstream  Mach 
number  and  decreases  with  an  increase  of  Mach  number,  The 
detailed  quantities,  such  as  location  of  the  RSP,  velocity  on 
dividing  streamline,  etc.,  arc  only  weak  functions  of  Mach 
number. 

4)  The  base  pressure  is  a weak  function  of  the  upstream 
boundary-layer  thickness.  All  oilier  quantities  such  as  cen- 
terline Mach  number,  etc.  are  strong  functions  of  Hence, 
measurement  of  the  latter  quantities  can  shed  more  light  on 
and  help  in  improving  the  corner  model. 

5)  At  low  base  bleed  rates,  the  theory  shows  a much  smaller 
base  pressure  rise  compared  with  experiments.  The  inaccuracy 
of  the  boundary-layer  equation,  in  representing  the  region 
close  to  the  base,  may  be  the  major  reason  for  this  result. 

In  conclusion,  the  method  developed  here  is  com- 
putationally fast,  provides  adequate  details  of  the  near  wake, 
and  hence  should  be  quite  useful  for  preliminary  design 
purposes. 
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Abstract 

A parametric,  study  for  an  external  burning  propulsion  system  was  carried 
out  to  investigate  potential  performance  and  application  of  such  a system. 

Ihe  study  used  a simplified  analysis  for  an  external  burning  zone  incorporated 
in  a realistic  near  wake  axisymmetric  supersonic  turbulent  base  flow  model 
based  on  the  Crocco  and  Lees  approach.  Ihe  axisymmetric  base  flow  model  with- 
out combustion  has  previously  been  shown  to  agree  well  with  experiment.  The  com- 
bustion zone  is  treated  as  a quasi  one -dimensional  zone,  and  is  considered  both  shear 
free  and  non-conducting.  With  these  major  assumptions,  it  is  shown  that  it  is  possible 
to  achieve  very  high  base  pressures  using  external  burning.  The  specific  impulse 
values  obtained  are  lower  than  reported  for  the  two-dimensional  case.  In  many 
cases  the  performance  is  less  than  usually  achievable  by  conventional  rocket 
propulsion,  if  current  fuel  rich  gas  generators  provide  the  fuel  for  the  ex- 
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ternal  burning  propulsion  system.  Recent  experiments  carried  out  using  com- 
pression surfaces  to  simulate  external  burning  confirm  the  order  of  magnitude 
of  specific  impulse  obtained  by  the  theory. 
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•r  turbulent  shear  stress 


rl’T2  corrections  due  to  non-zero  velocity  gradients  caused  by  entropy 
gradients  at  the  edge,  of  the  wake 

X heat  distribution  parameter 


Subscripts 


value  at  the  base 


Bsi- 


value  in  external  flow 


incompressible  quantities 

value  at  centerline  ; stagnation  quantity 

value  at  large  distance 

upstream  condition 


Superscript 


average  quantities  in  one-dimensional  zone 


1.  Introduction 

There  are  many  flight  missions  which  require  or  may  benefit  from  substan- 
tial drag  reduction  or  elimination  during  a portion  of  a supersonic  trajectory. 
It  has  been  recognized  for  some,  time  that  external  burning  during  supersonic 
flight  can  produce  compression  waves  which  can  be  translated  into  propulsive 
forces  on  a flight  vehicle.  This  external  burning  concept  is  different  from 
the  relatively  better  explored  base  burning  concept In  the  externaL  burning 
method,  the  combustion  mostly  takes  place  in  the  inviscid  flow  adjacent  to  the 
viscous  wake,  while  in  the  base  burning  method,  the  low  velocity  combustible 
mixture,  injected  through  the  base,  burns  in  the  viscous  wake.  There  is  a 
limitation  on  the  maximum  base  pressure  rise  achievable  in  the  base  burning 
concept.  There  is  no  such  limitation  in  the  external  burning  concept.  Hie  ex- 
ternal burning  method  has  also  been  claimed  to  be  simple  in  design  and  to  give 
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high  performance. 

The  external  burning  concept  has  initially  been  investigated  in  the 

(2) 

two  dimensional  case  by  Strahle,  using  the  original  Crocco  and  Lees 
theory of  the  near  wake  flow.  He  showed,  (i)  base  pressures  much  higher 
than  the  ambient  pressure,  and  (ii)  reasonable  specific  impulse  values  can 
be  achieved  using  this  concept.  Smithey^  and  recently  Neale  et  al.^, 
carried  out  experiments  with  axisymmetric  bodies,  using  compression  surfaces 
to  simulate  external  burning,  at  Mach  2 and  3,  respectively.  Their  results 
confirmed  the  first  finding  of  high  attainable  base  pressure.  But  much  smaller 

(4) 

specific  impulses  were  indicated  than  anticipated,  Smithey  also  applied  the 
external  burning  concept  to  the  axisymmetric  case,  using,  again,  the  original 
Crocco  and  Lees  theory.  This  original  theory  is  known  to  contain  an  inadequate 
prediction  of  the  length  scales  of  the  near  wake.  Recently,  the  more  detailed 
theory  of  Alber  and  Lees^6\  applicable  to  the  two-dimensional  turbulent  case, 
was  extended  by  the  present  authors  (and  also  independently  by  Peters  etall8^) 
to  the  axisymmetric  case;  it  was  shown  to  accurately  predict  the  neor  wake  de- 
tails. In  the  present  paper,  this  theory  has  been  modified  to  treat  external 
burning  propulsion.  Calculations  are  carried  out  to  predict  the.  performance 
of  such  a system  more  realistically. 

A short  outline  of  the  theory  without  combustion  and  the  basic  flow  model 
used  in  this  theory  arc  presented  in  Section  II.  For  more  details,  one  is  re- 
ferred to  the  back  up  paper  of  Reference  (7)  or  to  Reference  (9). 


Section  111  describes  the  details  of  the  heat  addition  zone  added  in  the 
inviscid  flow  region  of  the.  above  model.  Fuel  is  considered  fully  dispersed  in 
this  zone  and  there  is  no  mixing  across  bounding  streamlines.  The  entrainment 
of  the  air  before  combustion  can  be  taken  into  account  by  shifting  the  initial 
position  of  the  streamlines.  In  reality,  there  will  be  a diffusion  flame;  this 
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effects  of  a diffusion  flame.  On  either  side  of  the  combustion  zone  the 
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flow  field  is  treated  by  the  method  of  characteristics  (MOC),  Between  the 
near  wake  boundary  and  the  heat  addition  zone  the  MOC  is  rotational  to 
allow  for  injection  shock  entropy  rise  and  to  treat  the  shedding  of  the 
initial  boundary  layer.  Outside  the  heat  addition  zone  the  MOC  is  irrota- 
tional.  Hie  heat  addition  zone  is  displaced  outward  a sufficient  amount  so 
that  heated  vStreamlines  do  not  intersect  the  near  wake  boundary  before  the 
wake  closure  singularity  (Crocco-Lees  critical  point)  is  encountered  in  a 
downstream  numerical  integration.  This  allows  the  near  wake  to  be  treated 
under  an  adiabatic  assumption.  Thus,  the  analysis  of  the  inner  region  and 
the  corner  region  remain  unchanged,  except  for  slight  modifications  in  the 
program  to  take  into  account  the  possible  compression  instead  of  expansion 
at  the  corner. 

In  the  Appendix  some  computations , using  the  Mach  number  distribution 
provided  by  the  compression  surface  • of  the  experimental  set  up  of  Reference 
(5),  are  shown.  The  agreement  with  experiment  is  good,  although  the  base  pres- 
sures obtained  arc  about  15%  less  than  those  obtained  in  the  experiments. 

Section  IV  deals  with  the  results  and  discussion  of  the  present  external 
burning  near  wake  theory.  It  is  first  shown  that  it  is  possible  to  achieve 


high  base  pressure  rises  without  failure  of  the  present  computational  scheme. 
However,  the  total  temperature  rise  and  other  values  of  the  parameters  have 
been  limited  either  by  transonic  Mach  numbers  in  the  heat  addition  zone  or 
by  the  heat  addition  zone  mixing  with  the  viscous  wake.  Next,  a parametric 
study,  based  upon  the  parameters  obtained  through  dimensional  analysis,  is 
carried  out.  The  reference  conditions  chosen  for  carrying  out  this  study, 
although  arbitrary,  are  selected,  firstly,  to  give  a large  variation  of 
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various  parameters  before  any  of  the  above  mentioned  limitations  come  into 
the  picture.  Secondly,  they  are  chosen  to  avoid  large  gradients  in  field 
quantities;  large  gradients  require  finer  step  size  and  hence  more  computer 
time.  With  this  study,  the  importance  of  the  various  parameters  are  delineated. 
Other  than  the  upstream  Mach  number  and  the  initial  boundary  layer  thick- 
ness, which  are  more  or  less  design  inputs,  it  is  essenti" 1 to  have  an  input 
as  to  how  to  distribute  the  combustion  to  obtain  the  most  efficient  utiliza- 
tion, and  this  also  has  been  discussed  in  a general  manner  in  this  section. 


II.  Highlights  of  the  theory  of  the  Near  Wake 
Without  Combustion 

This  theory,  which  has  been  presenter  in  Reference  (7),  is  an  extension 
of  the  work  of  Alber  and  Lecs^  on  two-dimensional  turbulent  supersonic  near 
wakes.  The  major  points  of  variance  arc  the  treatment  of  the  comer  region, 
the  detailed  form  of  the  eddy  viscosity  model,  and,  of  course,  the  treatment 
of  the  external  flow.  The  model  used  to  represent  the  near  wake  flow  field 
is  as  shown  in  Tig.  1.  A uniform  supersonic  flow  at  constant  pressure  and 
Mach  number,  with  a finite  boundary  Layer  thickness,  approaches  the  base  of 
an  axisymmetric  body.  It  undergoes  an  expansion  near  the  corner.  Hie  viscous 
part  beyond  the  base  consists  of  (1)  the  recirculatory  region  and  (2)  the 
shear  layer  region,  while  the  inviscid  part  consists  of  (3)  the  rotational 
layer,  which  merges  into  the  shear  layer,  and  (4)  the  irrotational  region. 

The  flow  diagram  of  the  approach  is  shown  in  Fig.  2.  The  input  variables 
are  Mej,  6j/R  and  I.  Hie  unknown  of  the  problem  is  base  pressure,  and  to  start 
with,  this  is  assumed.  Hie  complete  flow  field  is  assumed  to  be  adiabatic  with 
turbulent  Prandtl  number  equal  to  1 and  to  obey  the  perfect  gas  laws.  For  the 
sake  of  analysis,  the  flow  field  is  divided  into  two  parts,  viz. , the  corner 


flow  region  and  Clio  flow  downstream  of  the  base.  The  solution  of  the  corner 
region  provides  the  Initial  conditions  for  the  flow  downstream  of  the  base. 
Specifically,  the  thickness  of  the  boundary  layer  in  which  the  viscous  forces 
remain  predominant  is  determined  hy  this  solution,  ibis  region  is  solved  by 
an  approximate  model.  Ihe  major  stipulations  in  this  model  arc  that  the  in- 
itial boundary  layer  can  bo  represented  by  a l/7ch  power  profile,  — = 

, and  there  i.s  an  isentrooic  expansion  in  the  outer  stveaintubes, 

of/o^ 

Downstream  of  tbe  base,  the  outer  inviscid  regime,  consisting  of  regions 
(3)  and  (4),  is  solved  vising  an  approximate  method  of  characteristics  (AMC), 
Ibis  method  is  a modification  of  Wcbb's^^  method,  and  is  consistent  with 
the  accuracy  and  spirit  of  this  approach,  The  inner  viscous  regime,  consist- 
ing of  region  (1)  and  (2),  is  represented  by  integral  boundary  layer  equations. 
The  governing  equations  used  to  represent  this  region  in  the  present  formula- 
tion are  given  below. 


Mass  conservation 


d . . d6 
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Momentum  conservation 
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Mechanical  energy  conservation 
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In  addition  to  these,  one  more  equation  Is  needed  near  the  base  iior  the 
chosen  velocity  profiles.  This  is  taken,  for  simplicity,  as  the  centerline 
moinen turn  e <|uu t ion . 
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As  the  flow  is  turbulent:,  some  modelling  of  the  shear  stress  is  re- 
quired, Here  shear  stress  is  evaluated  using  the  turbulent  viscosity  con- 
cent, the  details  of  which  are  given  in  Reference  (7) . The  integrals  in  the 
above  equations  are  obtained  in  closed  form  using  a compressible-incompressible 
transformation.  The  solution  is  started  from  the  base  using  two  parameter 
Green's  profiles  which  give  a good  representation  of  the  velocity  profile 
near  the  base,  with  or  without  uniform  base  bleed.  Those  are  given  by 
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The  use  of  these  profiles  results  in  six  unknowns,  viz,,  p (or  u^),  A, 
v , v,  lu/6^  and  S . The  oxtornuL  flow  solution  provides  v^(or  Q^)  and  S^, 
while  the  solution  of  Eqs.  (!)  - (4)  provides  the  remaining  unknowns.  Away 
from  the  base,  the  velocity  profiles  are  represented  by  the  Kebota-Recvcs- 
Buss^-^  one  parameter  profiles.  As  the  number  of  unknowns  arc  reduced  by 
one  using  these  profiles,  only  Eqs.  (1)  - (3)  ore  sufficient  to  represent 
the  inner  region.  The  uniqueness  of  the  solution  in  this  approach  comes  from 
the  occurrence  of  the  Crocco-bees  singular  point  downstream  of  the  wake.  The 
physical  explanation  of  this  singularity  is  that  the  near  wake,  which  is  sub- 
sonic in  the  mean,  becomes  supersonic  in  the  mean  after  this  point.  Mathe- 
matical ly,  this  singularity  imposes  an  extra  constraint  on  the  problem.  And 
the  base  flow  problem  turns  out  to  be  a boundary  value  problem,  with  base 
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pressure  as  die  eigenvalue,  and  u unique  solution  of  the  near  wake  la  ob- 
t:a  i.ned  , 

T!hs  strong  point.;',  of  this  theory  are  the  small  computer  time  required 
to  solve  the  problem  and  accurate,  prediction  of  tho  near  wake  details. 


ill.  Theory  with  External  Burning 

Tito  typo  of  external  burning  considered  here  is  the  concept  which  allows 
compression  waves  to  impinge  upon  the  near  wake  behind  a bluff-base  body  and 
to  thereby  raise  the  buse  pressure.  In  order  to  gain  a feel  for  the  perform- 
ance and  to  bring  out  tho  essential  features  of  an  external  burning  propul- 
sive system,  a simple  model  is  used  to  represent  the  external  combustion.  The 
external  burning  is  limited  to  the  inviscid  flow  adjacent  to  the  near  wake 
ns  shown  in  Fig,  3,  Fuel  is  considered  fully  dispersed  in  this  zone  and  there 
is  no  mixing  across  boundary  streamlines.  The  heat  addition  zone,  is  displaced 
outward  u sufficient  amount  so  that  heated  streamlines  do  not  intersect  the 
near  wake  boundary  before  the  wake  closure  singularity  (i.voceo-i.ee.s  critical 
point)  is  encountered  in  the  downstream  numerical  integration.  This  allows 
the  near  wake  to  bo  treated  under  an  adiabatic,  condition,  just  like  in  the 
non-burning  case.  On  either  side  of  the  heat  addition  zone  the  flow  field  is 
treated  by  the  AMC  as  in  Section  1.1,  between  the  near  wake  boundary  and  tho 
heat  addition  zone  the  AMC  is  rotational  to  allow  for  injection  shock  entropy 
rise  and  to  take  care  of  the  shedding  of  tho  initial  boundary  layer,  Tho  heat 
addition  zone  itself  is  treated  l>y  one -dimens iona 1 gasdynamics,  but  with  a 
centrifugal  correction.  The  governing  equations  used  to  describe  this  zone 
are  the  one -dimensional  energy  equation 


'"l-D  Af0  c 


(f>) 
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the  one-dimensional  continuity  equation 


mi  _r»  p A V 


the  strearowise  momen turn  equation 


- f = *1-D  d" 


and  the  momentum  equation  normal  to  the  streamlines 


o u2  is it 


ISq.  (6),  upon  manipulation,  can  be  rewritten  as 


To  d* 


T W 

L_  / Vt2_  ^ _JL 

1 + F \ T0  / Cp'ro,o 


This  gives  the  distribution  of  total  temperature  knowing  the  heat  release 


rate  function.  However,  in  the  present  calculations,  the  distribution  of 


total  temperature  is  assumed;  the  overall  energy  equation  is  used  to  obtain 


the  fuel-air  ratio. 


cJ’/iif 

pot 


1 -c  AT  /Hf 

pot 


(10a) 


Eqs.  (7),  (8)  and  (9)  upon  differentiation  and  simplification  give 

2 . 2 (ru  Ifu  *ll  ) 

1 dp  l+(y-l)M  1.  dM  dX  4 dx 

■“  % ^ ' in—  i — ■ r\  a 


5 dx  lt^  M2  « 


1 dp 


p dx 


i »w  ^2  dT 

1 dM  y M o _ n 

- J-  + " J"  " 


-2  M dX  2 To  A* 


srrt-J i : V > 


j_  ■-»  2i  1 

YH2  5 V't 


In  addition  to  these,  the  geometrical  relations  needed  are 


r - r , 
u u-i 


r-tWVl 


= tan  0 


tan  0, 


Figure  4 shows  the  points  referred  to  in  the  wake  flow  in  the  above  equations 
and  aiso  the  points  referred  to  later,  p and  0 are  defined  as 

pu+>'t 


e - S^t 


Also,  as  there  cannot  be  a pressure  difference  or  flow  angle  difference 

across  a streamline,  P^=PMi2  * Pu  = P1  ' ' = ^12  and  6u  = Q1 ' * However> 

’■locity  and  temperature  differences  across  streamlines  are  allowed;  hence 

uu. „ may  not  be  equal  to  u„,  etc. 

Ml  i.  £ 

The  unknowns  involved  here  are  p , 0 , M , p , 0 , M,,r. ,r  and  x . Eqs,(12)- 

U U U 1 u 

(15),  with  the  relations  provided  by  the  method  of  characteristics  are  suffici- 
ent to  solve  for  these  variables.  The  solution  scheme  used  is  a single  vari- 
bic  (0j'  ) iteration  scheme,  since  it  gives  reasonable  accuracy,  and  is  ex- 
plained below: 

(i)  Locate  point  Ml 2 by  extrapolation  of  ou  and  0,  and  using  the  mass 
conservation  equation  up  to  the  rotational  layer. 
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(ii)  Extrapolate  v^2  ^rom  the  known  values  of  v at  points  Ml  and  Ml-1. 
Calculate  ^^2*  auc*  t^en  ^Mi2  us^n8  isentropic  relations. 

(iii)  Obtain  Gj^  usin8  the  relation  valid  for  left  running  characteristics 
in  the  rotational  flow,  i.e., 


9M12  “ VM12  VM1  + ^1 


1 Srn  6M1  ^*0 Mi 2. Ml 
^11  rMl 


vm 


SM12  “ SM1 


VR  oHU  Mm1' 


(iv)  Obtain  dQ/dx  using  Eq.  (9),  and  the  new  G.  From  Eq.  (12),  find  Qu.  Also 
using  geometrical  relations  (10)  and  (11),  calculate  r^. 

(v)  Find  dM/dx  using  Eqs.  (7)  and  (8),  viz. 


dM 

dx. 


- M 


1 + M2 
1 - M2 


2 I' 


2 2 
dr  z - dr/ 
u £ 


dx 


dx 


2 2 
U ru  “ rt 

9 dr 

( 1 + v m ) 


dx 


(17) 


i 


Then  obtain  dp/dx  using  F.q,  (8),  and  hence,  new  values  of  p and  p^  . 

(vi)  Next  obtain  M/using  the  isentropic  relation, and  calculate  . 

(vii)  Obtain  q/  using  the  relation  along  right  running  characteristics  in 
the  rotational  region,  i.e.. 


Sin  0j / ^ 

ei ' = vi  '-i + ei  '-1  - vi ' + 5in  h •- 1 


A V , i'-i 


(viii)  If  the  above  value  of  0^ / is  not  the  same  as  that  of  0^  obtained  in 
step  (iv),  guess  a new  value  of  p^^,  ancl  °btaine  the  solution  by 
iteration  by  repeating  steps  (iii)  to  (vii). 
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Tlvc  specific  impulse  is  obtained  using  its  definition 


which  gives 


«P 


Decrease  in  drag 

Mass  flow  rate  of  fuel 


^,pb  Ab 
= » 1 •■■  1 • 

m£  g 
f o 


^b/pl  P1  Ab 

LsP  = F/<l+F)  mi_D  go 


(18) 


Finally,  two  major  points  of  interest  should  be  noted  in  the  present 
analysis.  First,  consider  Eq.  (17).  The  second  term  in  the  parenthesis 

mostly  dominates  at  least  by  one  order  of  magnitude  the  contribution  of 

- dM 

the  first  term  for  large  M.  But  as  more  and  more  heat  is  added,  as  — 

is  negative  for  M > 1,  M moves  towards  1,  and  the  order  of  magnitude  of 
both  the  terms  becomes  equal.  Thus  the  approximation  of  one  dimensionality 
becomes  very  poor  near  M equal  to  1.  Secondly,  in  the  concept  of  external 
burning,  there  is  high  loss  of  energy  because  very  hot  gases  leave  the  sys- 
tem. 

IV,  Results  and  Discussion 

The  present  calculations  are  made  with  a fuel  of  comparatively  low  calo- 
rific value,  10,000  Btu/lb,  as  this  value  is  easily  achievable  from  current 

(13) 

fuel  rich  gas  generators  . The  freestream  temperature  is  assumed  to  be 
519°R.  The  combustion  Is  assumed  to  be  1007.  efficient.  The  maximum  ratio  of 
fuel  to  air  that  is  permitted  is  the  stoichiometric  fuel-air  ratio,  since 
no  mixing  or  diffusion  flame  is  allowed  in  the  present  model.  However,  a 
rough  computation,  using  a simple  fuel,  showed  that  this  limit  is  reached 
after  the  M = 1 point,  where  the  present  computations  are  invalid.  The 
molecular  weight  etc.  of  the  mixture  is  taken  as  that  of  the  air. 

In  order  to  check  whether  or  not  the  inclusion  of  the  one-dimensional 
zone  would  substantially  degrade  the  original  theory  of  the  near  wake,  a 
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tow  no-heat  addition  cases  with  various  dimensions  of  the  one-dimensional 


zone  were  computed  and  compared  with  the  original  method  using  the  AMC  in 
the  outer  flow.  The  agreement  between  the  centerline  pressure  distributions 
of  the  approximate  and  exact  solutions  is  within  57>  everywhere,  showing  that 
a one-dimensional  treatment  of  the  heat  addition  zone  should  not  yield  un- 
acceptable error.  However,  this  variation  must  be  kept  in  mind  when  inter- 
preting some  of  the  results  showing  the  effect  of  various  parameters. 

Figures  5 and  6 arc  the  results  of  the  computations  carried  out  close  to 


the  M = 1 limit,  with  total  temperature  rise  (or  fuel-air  ratio,  see  Eq.(lO-a)) 


as  the  parameter.  These  runs  confirm  the  finding  of  previous  research  worker^*^’^ 


that  high  base  pressure  rise  is  possible  using  this  concept.  Also,  this  figure 
gives  a feel  for  the  specific  impulse  values  that  can  be  expected  at  practical 
base  pressure  rises  with  this  method.  The  specific  impulse  decreases,  i.e., 
the  system  becomes  less  efficient, both  with  increase  of  the  Mach  number  and 
increase  of  the  fuel-air  ratio.  The  minimum  value  of  the  specific  impulse 


^ T 

shown  in  Figure  5 is  about  86  secs,  occurring  at  Mach  3.0  with  ^ o 


To,c 


- 2.27, 


with  a corresponding  base  pressure  value  equal  to  1.85  times  the  ambient 
pressure.  The  maximum  value  of  I is  about  198  secs,  occurring  at  Mach  2.0 

with  a corresponding  base  pressure  value  equal  to  0.739 


with  ^To 


= 0.267, 


O,® 


times  the  ambient  pressure.  The  rear  stagnation  pressure  and  the  base 


pressure  rise  with  the  increase  of  fuel-air  ratio,  but  there  is  not  a one-to- 
one  correspondence  in  their  rises.  For  a fixed  Mach  number,  this  is  mainLy 
due  to  the  movement  of  the  rear  stagnation  point  (RSP).  The  location  of  the 
RSP  tends  to  reach  a limiting  value  at  higher  fuel-air  ratio,  and  then 


fip  /Ap,  becomes  almost  constant.  This  ratio  is  also  affected  by  Mach  number, 

Iwi  D 


and  other  parameters  which  move  the  location  of  the  RSP. 
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Next  a parametric  study  was  conducted  to  delineate  the  importance  of 
various  design  variables.  From  dimensional  analysis,  the  specific  impulse 


can  be  written  in  terms  of  the  non-dimensional  parameters  as  follows 
•^sp&o  r„  xi,  1*1  AT, 


f (v  M -1  rii  ^ XiL  LL  .^52_  } 

c , *'X»  el’  R ’ R ’ R * R ’ R ’ T ’ H , ' 

'-el  o.®  f 


To  conduct  this  parametric  study,  the  reference  or  design  condition  is 
chosen  on  the  basis  that,  firstly,  it  should  allow  a study  of  a large  vari- 
ation of  parameters,  and  secondly,  the  computer  time  for  this  study  should 
be  low.  For  a high  heat  addition  and  a high  Mach  number  design  point,  the 
gradients  occurring  in  the  flow  field  are  large,  and  a smaller  step  size 
and  hence  more  computer  time  is  required.  Hie  present  study  is  made  at  medium 
heat  addition,  uniform  combustion,  with  the  following  values  of  the  parameters: 


Mcl  ' 2*° 

Xy/R  = 1.5 

Oy  Xj,)/R  = 2. a 

rL/R  = 1.5 

ru/R  * 1.9 

6j[/R  = 0.181 

AT  /T  = 0.398 

O O ,oo 

c T /H,  = 0.01253 
p <»  f 

Table  1 shows  the  effect  of  the  heat  distribution  parameter.  Three  types 
of  heat  distributions  were  studied,  and  these  correspond  to  (a)  a combustion 


zone  which  decreases  in  strength  with  distance,  i.e,. 


dT 


= K.(l- 


VX1l 


>. 


dx  1 

(b)  a constant  combustion  rate  throughout  the  heat  addition  zone,  i.e,, 
dT 

— — = K„  , and  (c)  a more  realistic  type  of  combustion  - a slow  rise  in  the 
dx  1 
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Table  1.  Effect;  of  heat  distribution  parameter 


Representation 


I M ..  g 
sp  el  o 


intensity  of  combustion  representing  the  ignition  stage,  followed  by  uniform 
combustion  and  then  a tapering  off  stage.  As  shown  in  Table  1,  there  is  virtu- 
ally no  effect  of  the  distribution.  However,  it  should  be  noted  that  this  inter- 

pretation is  correct  only  if  the  combustion  zone  is  small,  so  that  all  the  com- 
pression waves  hit  in  the  initial  portion  of  the  near  wake.  If  the  combustion 
zone  has  a large  length,  type  (a)  will  be  the  best  since  there  will  be  less 
wasted  heat. 

The  effect  of  freestroam  Mach  number  on  specific  impulse  and  base  pressure 

rise  is  show  in  Fig.  7.  I M , g /c  , and  Ap,  /(p,  M ,)  arc  plotted  instead  of 

sp  el  o cl  rb  1 el  r 

I g /c  , and  Ap,  /v-,  since  they  are  less  dependent  on  the  Mach  number.  I 
sp  o el  rb  1 sp 

falLs  from  252  secs,  to  110  secs,  as  M , is  increased  from  1.5  to  3,5.  The 

el 

variation  of  some  typical  features  of  the  near  wake,  xRSp  and  PRgp>  are  also 
plotted  on  the  same  figure.  The  variation  of  these  is  ulmost  opposite  to  what 
is  obtained  without  heat  addition. 

Figure.  8 shows  the  variation  of  various  variables  with  upstream  boundary 
layer  thickness.  The  effect  of  the  heat  addition,  in  general,  is  to  nullify  the 
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effect  of  boundary  layer  thickness  on  base  pressure  without  combustion,  The  specific 
impulse  reaches  an  asymptotic  value  after  a moderate  thickness  of  the  boundary  layer. 
The  effect  of  variation  of  the  combustible  mass  by  varying  the  upper 
stream  line  position  is  shown  in  fig.  9,  The  value  of  1 increases  from 

302  to  142  sees,,  and  n,  increases  from  0.638  to  0.9313  as  r /R  (.<;  raised 

. b . H 

lUi  jv  \\\\  tx 

from  1.6  ( ” — — « 0,31)  to  2.4  ( ~ — t = 3.51).  AlWp/Ap.,,  the 

Pel  uclAb  °cl  cl  Ab  KbP  b 

shape  parameter  of  the  centerline  pressure  curve,  almost  remains  constant 
(not  shown),  while  X.ggp/R  increases  with  an  increase  of  r^/R. 

Figure  10  shows  the  variation  of  specific  impulse  parameter  and  pressure 
rise  parameter  with  axial  location  of  the  heat  addition  zone.  The  flat  maxi- 
ma in  these  curves  verify  the  results  obtained  with  the  variation  of  heat 
distribution  parameter,  x»  However,  as  the  heat  addition  zone  is  moved  far- 
ther and  farther  away  from  the  base,  part:  of  the  compression  due  to  external 
burning  affects  the  flow  downstream  of  the  wake  critical  point,  and  is  in- 
effective in  a base  pressure  alteration.  Finally,  a condition  is  reached  when 
the  presence  of  external  burning  is  not  Celt  at  all  by  the  near  woke,  and  the 
same  near  wake  details,  as  without  combustion,  are  obtained.  This  figure  also 
shows  the  variation  of  xRSp/R  and  Ppgp/Pp  xg$p  varies  linearly  with  x^  in  most 
of  the  portion,  as  expected  from  the  experiments  with  compression  surfaces. 

Figure  11  shows  that  there  is  an  optimum  rate  of  burning  the  fuel  (or  heat: 
addition  rate).  The  maximum  is  flat,  as  would  be  expected  from  the  studies  of 
variation  of  x ;nid  x^,  The  sharp  fall  in  performance  parameter  after  some 
stage  is  again  due  to  the  fact  that  part  of  the  heat  addition  becomes  useless. 
The  characteristic  properties  of  the  centerline  pressure  distribution,  x^^/R 
and  Al>j,gp/Al>|} , have  an  almost  similar  shape. 

The  effect  of  the  radial  location  of  the  heat  addition  zone  on  the  perform- 
ance, as  shown  in  Fig.  12,  is  very  weak.  This  is  again  due  to  the  basic  mechan- 
ism of  raising  the  base  pressure  in  an  external  burning  system.  The  compression 
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Figure  10,  Effect  of  Axial  Location  of  Combustion  Zone* 
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waves  caused  by  the  heat  addition  meet  the  viscous  wake  before  the  critical 
point  for  almost  the  whole  range  of  r^/R  shown  in  the  figure,  k^gp  varies 

linearly  with  r^  (as  would  be  expected  from  the  study  of  variation  of  x^) , 

while  PRgp  remains  almost  constant. 

Figure  13  shows  the  effect  of  ^o  , ^ . Curves  very  similar 

T ^ ~ t T ' * 

o,®  p O,co 

to  those  presenting  the  effect  of  the  change  of  combustible  mass  are  obtained. 

I varies  from  250  to  189  secs,  and  p.  /p.  varies  from  0.613  to  1.011  as 
sp  D L 

AT  /!  changes  from  0.07  to  0.94.  The  location  of  the  RSP  and  ApHQp/Ap. 

O O Ko 1 b 

approach  asymptotic  values  with  high  values  of  ATq/To  a). 

Finally,  the  variation  in  performance  with  the  freestream  temperature  is 

shown  in  Fig.  14.  This  variation,  however,  can  be  easily  predicted.  The  base 

pressure  is  independent  of  T or  Reynolds  number  when  the  flow  is  turbulent 

00 

(except  from  the  weak  effect  which  comes  through  the  boundary  layer  thickness). 

Using  this  information  with  the  definition  of  I as  given  by  Eq.  (18),  it  can 

sp 

be  shown  that  1 “ ; • 

From  the  above  study,  the  question  of  how  small  variations  (intentional 
or  unintentional)  from  the  reference  point  will  effect  the  performance 
can  be  answered.  Some  parameters  like  x > x^ , (x^-x^)  ant*  r^  have  previously 
been  shown  to  have  secondary  effects  on  the  performance  for  reasonable  design 
conditions.  Figure  1 5 is  the  plot  of  the  specific  impulse  parameter  against  the 
total  energy  parameter,  and  is  obtained  from  Figs.  9 and  13.  It  shows  that 
for  the  previous  reference  conditions,  it  is  more  efficient  to  change  the 
fuel-air  ratio  if  a little  higher  base  drag  reduction  is  required,  and  to 


change  r^/R  if  a little  lesser  base  drag  reduction  is  required  for  the  same 
projectile,  same  Mach  number  and  altitude  of  flight.  However,  no  optimiza- 
tion studies  have  been  carried  out. 


V.  Conclusions  and  Comments 

The  present  study  shows  that  the  performance  of  an  external  propulsion 
system  is  strongly  dependent  only  on  the  upstream  Mach  number,  fuel-air  ratio 
the  combustible  mass  and  the  fuel  calorific  value.  Also  the  system  becomes 
more  efficient  with  a cold  freestream.  For  preliminary  design  of  this  system, 
reasonable  values  of  heat  distribution  parameter,  axial  and  radial  locations, 
and  length  of  the  heat  addition  zone  can  be  assumed,  as  their  effect  on  the 
performance  is  only  secondary.  Thus  for  a given  projectile  flying  at  a given 
Mach  number  and  altitude,  and  given  base  pressure  rise,  the  optimization 
problem  is  reduced  to  the  selection  of  r^/ll  and  ATQ / Tq  oj.  One  rule  of  thumb 
which  came  out  of  these  studies,  and  should  be  helpful  in  extrapolating  ex- 
perimental results,  is  that,  except  for  very  low  heat  additions,  base  pres- 
sure remains  constant  with  Mach  number  awl  upstream  boundary  Layer  thickness 
for  the  same  values  of  other  parameters. 

live  question  whether  the  present  values  indicate  an  upper  limit  of  per- 
formance is  unanswered.  To  answer  this  question,  one  has  to  consider  two 
effects,  viz. , the  losses  due  to  viscosity  and  the  mass  entrainment  in  the 
heat  addition  zone  due  to  mixing.  An  improvement  in  modelling  of  the  external 
burning  zone  is  required  to  answer  this  question  satisfactorily. 

Finally,  the  present  calculations  show  that  although  high  net  thrust  can 
be  obtained  using  the  external  burning  method,  this  method  is  not  efficient 
unless  a high  Btu  fuel  is  used.  Its  performance  falls  short  of  that  of  con- 
ventional rocket  propulsion  systems.  However,  it  is  still  attractive  because 
of  the  simplicity  of  design.  Its  performance  probably  can  be  boosted  by  com- 
bining this  method  with  the  base  burning  method,  which  is  very  efficient  for 
low  base  pressure  rise  values.  More  experimental  and  theoretical  work  are 
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warranted  to  check  these  conclusions 
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Append  Lx 


(<j\ 

The  present  theory  was  compared  with  experiments  in  which  external 
compression  was  carried  out  by  modifying  the  shape  of  the  wind  tunnel  test 
section.  The  conditions  on  the  initial  Mach  wave  at  the  corner  were  obtained 
using  potential  fiow  theory,  ignoring  shocks  and  the  boundary  layer  thickness. 
The  effect  of  the  boundary  Layer  in  these  calculations  was  taken  into  account 
in  a similar  fashion  as  in  the  original  theory.  The  rotational  layer  was  ig- 
nored for  simplicity,  and,  also,  since  the  base  pressures  were  high  with  com- 
pression surfaces,  the  viscous  forces  remain  predominant  in  whole  of  the  initial 
boundary  layer.  Figures  10  and  17  show  both  the  experimental  and  theoretical 
results.  The  agreement  between  theory  and  experiment  is  reasonable,  although 
the  base  pressure  obtained  is  about  1 57,  less  than  given  by  experiment.  Some  of 
the  difference  may  be  due  to  (i)  the  approximations  made  in  calculating  the 
conditions  on  the  initial  Mach  wave,  (ii)  the  slight  compression  caused  by 
the  boundary  layer  on  the  tunnel  wall,  and  (iii)  the  slightly  Improper  account- 
ing of  the  boundary  layer  thickness  effect. 
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Figure  16.  Computations  with  External  Compression: 

Effect  of  External.  Compression  Strength 
and  Locution  on  Ba.se  1'rcssvive. 
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